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摘要 : 几何 扭曲 (Geometric Distortion, GD) 是 电荷 耦合 器 件 ( Charge-Coupled Device, CCD) 
高 精度 测量 不 可 忽略 的 因素 ， 通 常 需 要 细致 的 观测 工作 导出 ， 这 必然 占用 宝贵 的 观测 时 间 。 
利用 丽江 观测 站 2.4m 望远镜 观测 密集 星 场 的 CCD 图像， 初步 发 现 了 图 像 几 何 扭 曲 随 天 顶 
距 (Zenith Distance, ZD) 的 变化 规律 。 有 具体 地 ， 本 文 参考 Gaia DR3 星 表 提供 的 位 置信 息 并 进 
行 数 据 归 算 ， 改 正 了 较 差 颜色 折射 (Differential Color Refraction, DCR) 效应 ， 求 解 出 不 同 观 
测 时 间 CCD 的 几何 扭曲 模型 ， 并 使 用 四 阶 多 项 式 拟 合 。 结 有 果 表 明 ，CCD 图 像 的 几何 扭曲 与 
观测 时 的 天 顶 距 存 在 一 定 的 相关 性 。 在 像素 和 方向， 几何 扭曲 模型 四 次 项 x ，x*y,，x y 和 
ay, SRR ax foxy 与 天 顶 距 存 在 较 明 显 的 线性 关系 ; By 方向 ， 与 天 顶 距 存在 线性 关系 的 
系数 项 为 四 次 项 x* 和 wy ”、 三 次 项 wy 和 二 次 项 xy 
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望远镜 光学 成 像 系统 的 几何 扭曲 会 导致 观测 的 CCD 图 像 中 目标 星 像 偏离 原本 的 像素 坐标 位 置 ， 

这 是 一 种 普遍 存在 的 现象 。 该 现象 不 仅 存 在 于 空间 望远镜 ， 也 存在 于 地 面 望 远 镜 ， 对 视 场 中 星 像 位 置 
( Whar EAN TY Bm ee], fn, CEMA es Epa ( Hubble Space Telescope, HST) 的 光学 成 像 系统 
~ 中 ， 视 场 存在 严重 的 几何 扭曲 ， 其 中 第 2 代 广 角 行星 照相 机 (Wide Field and Planetary Camera 2, 
= bis e geri 为 了 准确 求解 哈 勃 空间 望远镜 的 几何 扭曲 效应 ， 文 [1] 提 
出 了 通过 提高 线性 项 精度 来 改正 几何 扭曲 的 方法 ， 使 得 行星 相机 视 场 精确 到 0. 02 pixel, 4 (Wide 
l Field, ee 01 pixel。 在 地 面 望远镜 的 观测 中 ，2006 4E, X2 ftne eyes ESLA 
c 扭曲 求解 方法 应 用 于 欧洲 南方 天 文 台 2. 2 m 望远镜 ， 取 得 了 约 7 mas 的 位 置 测 量 精度 。 随 后 ， 文 [3] 
© 提出 了 一 种 基于 参考 星 表 位 置 测 量 信息 的 几何 扭曲 求解 方法 ， 提 高 了 地 面 望远镜 的 观测 准确 度 与 精 
度 ， 并 成 功 应 用 于 土 卫 九 (Phoebe) 的 位 置 测量 。 此 后 ， 该 方法 应 用 于 多 项 研究 ”,， 文 [5] 在 校正 几何 
扭曲 后 得 到 近 地 小 行星 Apophis 的 位 置信 息 。 在 Gaia 星 表 发 布 后 ， 该 方法 成 功 应 用 于 木 卫 六 
( Himalia)' 与 海 卫 一 (Triton) 的 位 置 测量 '"。 文 [8] 参 考 Gaia DR2 星 表 ， 并 从 准确 度 、 精 度 等 方面 系 
统 比 较 了 文 [1] 与 文 [3] 的 方法 。 随 着 最 新 的 Gaia DR3 参考 星 表 发 布 ， 文 [3] 所 述 方法 的 应 用 将 更 广 
泛 ， 能 够 更 精确 地 测量 天 体位 置 。 

文 L[1] 和 文 [3] 这 两 种 方法 在 求解 几何 扭曲 时 对 所 使 用 的 观测 资料 有 严格 要 求 ， 观 测 资料 需要 采 
用 不 同 取 向 或 指向 的 观测 模式 ， 均 勺 抖动 拍摄 ， 这 一 工作 通常 需要 充足 的 观测 时 间 和 一 定 的 观测 经 验 
才能 完成 。 然 而 ， 望 远 镜 的 观测 机 会 宝贵 ， 观 测 时 间 有 限 。 为 了 能 充分 利用 有 限 、 宝 贵 的 观测 时 间 , 
本 文 希望 从 天 项 距 方面 探究 云南 天 文 台 丽 江 观测 站 2. 4 m 望远镜 CCD 图 像 的 几何 扭曲 变化 规律 。 天 
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顶 距 是 天 体 与 天 顶 之 间 的 角 距 离 ， 天 顶 距 越 大 ， 大 气 折射 越 大 。 然 而 ， 大 气 折射 造成 天 体 的 视 方向 与 
真实 方向 产生 偏差 ， 影 响 天 体位 置 测量 ， 并 且 由 大 气 折射 造成 的 扭曲 项 是 光学 扭曲 的 重要 项 ， 文 [9] 
求解 了 大 气 折射 扭曲 项 对 暗 能 量 相机 天 体 测 量 定 标的 影响 。 因 此 ， 本 文 认为 天 顶 距 可 能 在 一 定 程度 上 
影响 几何 扭曲 模型 的 变化 。 本 文采 用 文 [3] 提 出 的 方法 求解 几何 扭曲 模型 ， 参 考 最 新 发 布 的 Gaia DR3 
星 表 进 行 数据 归 算 ， 并 考虑 较 差 颜 色 折 射 效 应 的 影响 ， 计 算出 不 同时 间 点 观测 的 几何 扭曲 模型 ， 试 图 
发 现 几何 扭曲 模型 随 观 测 天 项 距 变 化 的 规律 。 本 文 布 望 通过 该 规律 ， 能 够 为 后 续 的 几何 扭曲 求解 提供 
一 定 的 便利 ， 即 通过 较 少 次 数 定 标 场 的 观测 计算 CCD 图 像 在 其 他 时 间 的 几何 扭曲 ， 从 而 有 效 地 节约 
望远镜 观测 定 标 场 的 时 间 ， 降 低 望 远 镜 观测 定 标的 频率 。 


1 观测 资料 
R1 望远镜 和 CCD 参数 说 明 
本 文 采 用 云南 天 文 台 丽 江 观 测 站 24m 光学 Table 1 Specifications of the telescope and the CCD chip 


望远镜 在 2013 4p2 H4 H «6 Hosp mi E H Item Parameter 
NGC1664 与 近 地 小 行星 Apophis 观测 的 资料 探究 Focal length/cm 1920 
几何 扭曲 模型 的 变化 规律 ， 观 测 所 用 望远镜 和 F-ratio 8 
CCD 的 详细 参数 见 表 1。 本 文选 取 10 组 不 同 天 Diners panes aon 240 
顶 距 拉动 观测 的 密集 星 场 的 资料 ( 视 场 中 大 约 能 "MMC MEM Bus 
观测 到 700 颗 恒 星 ) ， 其 中 包括 66 帧 NGC1664 Size of pixel 13.5 um x13.5 pm 
星团 图 像 BRAT ARE ) pu VUE GE Size of CCD array ( effective) ] 900 x 1 900 
Apophis 图 像 ( 原 用 于 高 精度 位 置 测量 ) ， 具 体 的 
观测 资料 如 表 2 Scale/(" /pixel) - 0. 283 
表 2 观测 资料 概要 
Table 2 Observations overview 
ID_Set Obs date Target Frame/ Filter ZD (mean)/(°) Exp time/s 
1 2013-02-04 NGC1664 22/B 24. 7-28.7 (26.7) 40 
2 2013-02-05 Apophis 46/B 32. 3-35. 1 (33.4) 30 
3 2013-02-05 Apophis 47/B 32. 3-34. 1 (32.9) 30 
4 2013-02-05 Apophis 46/B 34. 0-38. 9 (36.2) 30 
5 2013-02-05 Apophis 45/B 38. 8-49. 5 (43. 7) 30 
6 2013-02-06 Apophis 22/B 31. 4-31. 8 (31. 6) 30 
7 2013-02-06 Apophis 22/B 31. 6-32.2 (31.9) 30 
8 2013-02-06 Apophis 45/B 32. 0-35. 2 (33.4) 30 
9 2013-02-06 NGC1664 22/B 46. 9-52. 4 (49.5) 40 
10 2013-02-06 NGC1664 22/B 44. 0-52.7 (47.6) 40 
2 Ji ik 


2.1 视 场 定 标 

对 于 每 幅 CCD 图 像 ， 采 用 如 下 步 又 进行 天 体 测量 归 算 。 首 先 ， 对 图 像 进行 平 场 与 本 底 校 正 ， 并 
裁剪 图 像 的 边缘 区 域 ， 裁 前 后 的 图 像 像素 大 小 为 1 900 x 1 900。 然 后 通过 二 维 高 斯 定 心 算法 测量 每 个 
星 像 的 像素 坐标 位 置 ， 并 将 像素 坐标 位 置 与 Gaia DR3 星 表 提 供 的 位 置信 息 进行 匹配 。 在 星 表 匹 配 过 
程 中 ， 为 避免 双星 系统 导致 的 位 置 误 差 ， 本 文 剔 除了 RUWE (Re-normalised Unit Weight Error) 值 大 于 
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2.0 的 恒星 ; 将 匹配 到 的 每 颗 恒 星 转换 到 观测 时 刻 的 丫 心 视 位置 ， 并 考虑 大 气 折 射 等 天 体 测 量 效应 ， 
再 通过 中 心 投影 公式 得 到 每 个 星 像 的 标准 坐标 ; 最后， 使 用 加 权 模 型 求解 四 阶 多 项 式 底片 常数 模型 ， 
得 到 每 幅 CCD 图 像 的 定 标 参数 ， 根 据 求解 的 底片 常数 模型 计算 每 个 星 像 的 观测 位 置 。 本 文 归 算 了 不 
同时 刻 的 观测 资料 ， 对 于 每 组 资料 ， 计 算得 到 了 每 颗 恒 星 的 残 差 (The Observed Minus Computed, O- 
C) 和 标准 差 。 
2.2 BEBE TTA RUT BCE 

由 于 光 在 空气 中 的 折射 率 与 波长 有 关 ， 不 同 波长 的 光 在 天 项 方向 受到 不 同 的 折射 ， 使 得 不 同 恒星 
的 光 穿 过 地 球 大 气 时 折射 为 不 同 的 光谱 。 因 此 ， 上 有 具有 不 同 光谱 的 天 体 成 像 位 置 受到 不 同 程度 的 大 气 折 
射影 响 ， 这 种 现象 称 为 较 差 颜 色 折 射 。 该 现象 可 能 导致 目标 天 体 的 位 置 测量 结果 包含 沿 天 顶 方向 的 系 
统 误差 。 本 文 使 用 的 观测 资料 均 采用 B 滤 光 片 拍摄 ， 该 滤 光 片 受 较 差 颜色 折射 影响 较为 严重 。 为 了 
求 出 更 准确 的 几何 扭曲 模型 ， 本 文 首 先 通 过 文 [10] 提 出 的 较 差 颜色 折射 改正 方法 消除 每 个 星 像 受 较 
差 颜 色 折射 的 影响 。 首 先 ， 参 考 Gaia DR3 星 表 恒星 的 位 置 与 测 光 信息 项 (BP-RP, 色 指 数 ) 匹配 图 像 
中 的 每 个 星 像 ， 根 据 星 表 中 的 天 体 测量 数据 ， 转 换 到 地 平 坐标 系 中 的 视 位 置 (包含 大 气 折 射 ) ， 得 到 
每 颗 测 量 星 的 地 平 高 度 和 方位 角 。 然 后 ， 通 过 加 权 的 四 阶 多 项 式 底片 常数 模型 归 算 得 到 每 颗 恒 星 的 残 
差 ， 使 用 一 阶 多 项 式 


DCR(a, b) =a + bcolortan( ZD) (1) 
拟 合 色 指 数 和 地 平 高 度 方向 的 残 差 之 间 的 关系 ,其 中 color 表示 参考 星 的 色 指 数 ; ZD 为 观测 资料 的 天 
TUE; a Ab 是 拟 合 参数 。 通 过 加 权 最 小 二 乘法 进行 拟 合 ,求解 较 差 颜 色 折射 改正 系数 a Mb, WE 
恒星 在 高 度 方 向 的 残 差 。 最 后 ， 对 每 个 星 像 的 像素 坐标 位 置 进行 校正 。 
本 文 对 10 组 观测 资料 均 进 行 了 较 差 颜色 折射 改正 。 图 1 给 出 了 较 差 颜色 折射 改正 前 后 资料 集 3 
在 天 顶 距 方向 的 残 差 随 色 指 数 的 变化 情况 ， 通 过 (1) 式 用 加 权 最 小 二 乘法 拟 合 ， 得 到 拟 合 参数 ， 蓝 色 
虚线 为 拟 合 结果 。 从 图 (b) 可 以 看 出 ， 改 正 较 差 颜色 折射 影响 后 ， 系 统 误 差 得 到 了 明显 改善 ， 表 明 较 
差 颜 色 折 射 对 天 顶 方向 的 影响 基本 改正 。 


— Fitted curve: a=30.8+0.2 mas, b=—48.4 0.2 mas/mag — Fitted curve: a—-0.3 30.2 mas, b-0.03-0.2 mas/mag 


O-C in ZD /(") 
O-C in ZD / (") 


0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 
BP-RP / mag BP-RP / mag 
(a) (b) 


图 1 (a) 较 差 颜色 折射 改正 前 资料 集 3 的 残 差 随 色 指数 的 变化 ， 蓝 色 实 线 表 示 (1) 式 的 拟 合 结果 ; 
(b) 较 差 颜色 折射 改正 后 的 结果 
Fig.1 (a) the (O-C) residuals change with the color index BP-RP for observation set 3 before DCR correction, 
the blue line represents the fitted result of Equation 1; (b) the result after DCR 


2.3 几何 扭曲 校正 
关于 几何 扭曲 求解 的 方法 , 文 [1] 通 过 自 定 标 方法 直接 将 图 像 的 像素 坐标 位 置 构建 为 自 洽 的 无 扭 
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曲 参考 框架 求解 几何 扭曲 ， 文 [3] 基 于 天 体 测 量 星 表 并 消除 相当 的 星 表 系 统 误 差 的 方法 求解 几何 扭 
曲 。 对 于 每 组 观测 资料 ， 通 过 本 文 第 2. 2 节 的 方法 改正 较 差 颜色 折射 影响 后 ， 通 过 如 下 步骤 求解 几何 
扭曲 。(1) 使 用 加 权 四 和 常数 底片 模型 进行 归 算 ， 归 算 步 又 见 本 文 第 2. 1 节 ， 计 算 每 幅 图 像 中 所 有 匹配 
星 像 的 残 差 值 (0-C) 。(2) 对 于 出 现在 WCV>3) 幅 不 同 的 CCD 图 像 中 的 星 像 ， 令 i 和 j 分 别 表示 其 中 
两 幅 不 同 的 图 像 ， 用 


dx, = Ax, - x X; 
e;cosD, e;cosD, (2) 
2 
A A 
dy, = Ay; - — Ay, + — dy, 
e e; 


oA BUE SIS i WMP 7 WR BR rp D ER ME B JL TRE SL HH (dx, dy; ) 和 (dx ，dy ) 之 间 的 数学 关系 ， 
其 中 , e = cose/p 为 四 常数 模型 求解 的 估计 参数 ; e Alp 分 别 为 CCD 的 取向 和 比例 尺 ; Ax 和 Ay 分 别 
为 星 像 的 测量 像素 位 置 (x,, y, ) 与 通过 使 用 估计 参数 的 四 常数 模型 间接 计算 得 到 的 该 星 像 的 像素 位 置 
(x, , y,) 2; D 为 天 球 切 平面 切 点 在 赤 纬 方向 的 赤道 坐标 。 通 过 求解 (2) 式 可 以 推导 该 星 像 像素 
> 位 置 的 几何 扭曲 。(3) 把 CCD 图 像 划 分 为 19 x 19 的 方形 区 域 ， 将 像素 坐标 位 于 同一 区 域 的 所 有 星 像 


e 的 几何 扭曲 求 平均 值 ， 并 把 该 平均 值 作 为 区 域 中 心 像素 坐标 位 置 处 的 几何 扭曲 。(4) 通过 双 线 性 插值 
© 计算 所 有 星 像 在 像素 位 置 处 的 几何 扭曲 值 ， 然 后 对 每 幅 图 像 所 有 星 像 的 像素 坐标 进行 校正 ， 得 到 初步 
© 改正 几何 扭曲 的 观测 数据 结果 。(5) 再 次 执行 上 述 求解 过 程 ， 得 到 几何 扭曲 增 量 ， 与 第 1 次 几何 扭曲 
一 模型 相 加 ， 得 到 第 2 次 几何 扭曲 模型 。(6) 重 复 迭 代 以 上 求解 过 程 ， 直 到 几何 扭曲 增 量 达到 预定 的 精 


i 度 要 求 。 
3 ”结果 分 析 


3.1 几何 扭曲 模型 变化 规律 
本 文 对 10 组 资料 进行 了 几何 扭曲 求解 ， 图 2 分 别 给 出 了 每 晚 典 型 的 几何 扭曲 矢量 图 (从 左 到 右 依 
次 为 资料 集 1，3，8) 。 我 们 使 用 四 次 多 项 式 
dx = 3 ax y 


0<itj<4,ijeN 


Q dy- > bx y 


Oxisj«4,i,je N 


(3) 


表示 几何 扭曲 模型 其中, x* 和 yy 为 输入 位 置 x 和 y 归 一 化 后 的 位 置 ， 范 围 在 [ -1, 1] 内， 
x = (x - 950)/950 
y = (y - 950)/950 ' 
由 归 一 化 位 置 可 以 更 容易 看 出 每 一 项 对 扭曲 值 影响 的 大 小 。 表 3 和 表 4 分 别 给 出 了 在 * 方向 和 Yy 方向 
不 同 天 顶 距 观测 资料 的 几何 扭曲 模型 多 项 式 的 四 次 项 、 三 次 项 和 二 次 项 的 系数 ， 表 的 第 1 行 依次 为 天 
顶 距 (以 (”) 为 单位 )、 四 次 项 系数 、 三 次 项 系数 和 二 次 项 系数 。 
图 3 和 图 4 分 别 给 出 了 在 * 方 向 和 Y 方 向 几何 扭曲 模型 的 四 次 项 、 三 次 项 和 二 次 项 系数 随 天 项 距 
的 变化 情况 ， 并 在 每 个 子 图 的 上 方 列 出 了 系数 项 、 系 数 的 均值 与 标准 差 。 我 们 初步 发 现 三 次 项 系数 的 
标准 差 比较 小 ， 系 数 比较 稳定 ， 而 四 次 项 和 二 次 项 的 某 些 系数 随 着 天 顶 距 的 变化 具有 一 定 的 规律 。 因 
此 ,我们 通过 线性 拟 合 分 析 各 项 系数 与 天 顶 距 的 关系 ， 图 3 和 图 4 中 蓝 色 实 线 为 线性 拟 合 的 结果 ， 表 
5 和 表 6 分 别 列 出 了 拟 合 结果 的 皮尔 森 (Pearson) 相关 系数 与 标准 差 (SD) 的 值 。 本 文 根 据 皮尔 森 相关 
性 分 析 各 项 系数 与 天 顶 距 的 相关 强度 ， 皮 尔 森 相关 系数 7 为 


(4) 
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图 2 几何 扭曲 模型 矢量 图 。 每 个 子 图 的 上 方 列 出 了 观测 日 期 以 及 以 像素 为 单位 的 扭曲 最 大 值 与 捏 曲 中 值 ; 每 个 子 图 
的 矢量 值 放 大 200 4& 
Fig.2 The vector graphs of GD models. At the top of each panel, the observational date, the maximum CD (Max) and the 


median GD (Med) are listed in units of pixels and a factor of 200 is used to exaggerate the magnitude of each GD vector 


X3 在 x 方向 几何 扭曲 模型 的 四 阶 多 项 式 系数 


Table 3 Fourth-order polynomial coefficients of geometric distortion models in the x-direction 


ZD/(°) Q40 a3, Ay i3 Q04 a39 Ay, Ay Q03 a% ay, Aon 
26.8 0.010 0.066 -0.069 -0.076 0.040 -0.936 -0.024 -0.907 0.006 -0.200 -0.439 -0.149 
33.4  -0.027 0.108 -0.002 -0.100 0.029 -0.918 -0.007 -0.869 -0.014 -0.238 -0.369 -0.184 
33.8  -0.023 0.076 -0.030 -0.118 0.026 -0.883 -0.005 -0.871 -0.017 -0.222 -0.355 -0.175 
36.3 0.007 0.120 -0.060 -0.118 0.001 -0.897 -0.028 -0.876 0.001 -0.221 -0.328 -0.143 
43.6 —0. 042 0.114 -0.027 -0.151 0.070 -0.899 -0.045 -0.887 -0.026 -0.189 -0.314 -0.206 
49.6 -0.189 0.173 -0.062 -0.169 -0.068 -0.881 -0.056 -0.864 0.000 -0.069 -0.241 -0.108 
47.7 —0. 077 0.198 -0.106 -0.103 0.050 -0.950 0.006 -0.901 -0.016 -0.241 -0.276 -0.159 
31.6 0.040 -0.051 0.023 -0.098 0.005 -0.867 -0.017 -0.849 0.001 -0.336 -0.382 -0.166 
31.9 0.049 0.057 -0.017 -0.065 0.016 -0.948 -0.026 -0.868 -0.020 -0.289 -0.432 -0.186 
33.4 0.031 -0.144 0.000 -0.136 -0.015 -0.879 -0.057 -0.893 0.000 -0.279 -0.319 -0.151 

表 4 在 ?方向 几何 扭曲 模型 的 四 阶 多 项 式 系数 
Table 4 Fourth-order polynomial coefficients of geometric distortion models in the y-direction 
ZD/(°) bao bs, by» bis bog bso by by bo, by b bo, 


26.8 
33.4 
33.8 0 
36.3 0 
43. 6 0 
49. 6 0.006 -0.018 -0.164 0 
T 0.082 -0.017 -0.073 0.079 -0.172 0.032 -0.916 -0.048 -0.931 -0.010 -0.193 -0.181 
6 0 
9 0 
4 0 


0.003 -0.033 -0.141 0.005 126 -0.016 -0.883 0.002 -0.934 0.032 -0.094 -0.242 


0.016 -0.002 -0.135 0.030 176 -0.008 -0.871 -0.015 -0.925 0.036 -0.138 -0.189 


0. 015 0.002  —0.107 117 -0.010 -0.871 -0.025 -0.896 0.058 -0.135 -0.262 


0. 001 0.015 -0.141 075 -0.037 -0.887 -0.004 -0.887 0.057 -0.121 -0.275 


0.012 -0.008 -0.110 102 -0.013 -0.879 -0.030 -0.938 0.050 -0.157 -0.278 


.045 -0.050 -0.838 -0.034 -0.903 0.071 -0.119 -0.296 
47. 


31. 
31. 
33. 


0. 067 0.029 -0.079 -0.073 0.045 -0.010 -0.851 -0.010 -0.837 0.108 -0.058 -0.354 


一 0. 007 0.008 -0.059 -0.025 0.179 -0.021 -0.883 -0.017 -0.821 0.022 -0.120 -0.234 


. 005 0.128 -0.003 -0.889 -0.011 -0.898 0.050 -0.142 -0.289 


0. 008 0.005 -0.136 
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图 3 在 x 方 向 几何 扭曲 模型 系数 随 天 顶 距 的 变化 情况 ， 蓝 色 实 线 为 线性 拟 合 得 到 的 结果 


Fig. 3 


blue lines are the results obtained from the linear fitting 


In the x-direction, the coefficients of geometric distortion models change with the zenith distance, 
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图 4 在 y 方 向 几何 扭曲 模型 系数 随 天 顶 距 的 变化 情况 ， 蓝 色 实 线 为 线性 拟 合 得 到 的 结果 
Fig.4 In the y-direction, the coefficients of geometric distortion models change with the zenith distance, 


blue lines are the results obtained from the linear fitting 
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Hop, n 为 样本 量 ; (X, YOM (X, Y) 分 别 为 两 个 变量 的 观测 值 和 均值 ;zx 和 ec; 为 样本 标准 差 。 相 
关系 数 绝对 值 大 于 0. 4 表示 两 个 变量 的 相关 程度 为 中 等 及 以 上 ， 本 文 取 相 关系 数 绝对 值 在 0.5 以 上 的 
为 关系 显著 。 在 x 方向,， KADEN ARONA x, xy, xy May, DRH Alay, TEYy JT, 
A& A sr E SCIO 2 VU x I xy? SERT cy MIT wy。 


表 5 线性 拟 合 结果 的 皮尔 森 相关 系数 表 6 线性 拟 合 结果 的 标准 差 
Table 5 The Pearson's correlation coefficient of Table 6 The standard deviation of the 
the linear fitting results linear fitting results 
Pearson's Pearson's Pearson's Pe: , 
Items ip S dn d Items pica S FS S Items SDinx SD iny || Items SDinx SD iny 
r nx rın y rinx rainy 
x* —0. 842 0. 555 xy  -0.173 0. 015 x* 0. 040 0. 032 xy 0. 022 0. 023 
x'y  -0.823 -0.026 xy 0.036  -0.858 xy 0. 030 0. 015 xy” 0. 019 0. 008 
xy  -0.510 —-0.088 y —0.288 -0.348 xy 0. 036 0. 037 y 0. 011 0. 040 
xy -0.685 0. 567 x 0.603 -0.165 xy? 0. 025 0. 035 x 0. 060 0. 033 
y —0. 185 0. 193 Xy 0.901 -0.587 y 0. 040 0. 053 Xy 0. 029 0. 031 
x 0. 032 0. 015 y 0. 269 0. 095 x) 0. 032 0. 023 y 0. 028 0. 054 
3.2 wW i£ 


从 系统 误差 角度 考虑 ， 儿 何 扭曲 现象 与 望远镜 的 重力 形变 和 大 气 折射 误差 有 关 。 地 面 望 远 镜 受 重 
影响 ， 望 远 镜 支 撑 结 构 、 光 学 元 器 件 发 生 弯 曲 或 变形 ， 导 致 地 面 观测 系统 的 扭曲 解 不 稳定 。 望 远 镜 
结构 的 重力 形变 产生 的 误差 通常 随 着 支撑 结构 的 仰 俯 角 变化 。 光 学 元 器 件 的 形变 一 定 程度 上 也 会 
引起 光学 系统 的 误差 。 地 面 望远镜 在 观测 目标 天 体 时 还 受到 大 气 折射 效应 影响 ， 从 而 导致 观测 到 的 目 
标高 度 比 实际 高 度 要 高 ， 这 种 效应 随 天 顶 距 的 增 大 而 逐渐 明显 。 虽 然 我 们 根据 观测 时 的 平均 温度 和 气 
压 考虑 了 大 气 折 射 模型 的 改正 ， 但 是 该 模型 的 改正 未 必 准 确 ， 例 如 湿度 、 风 力 、 观 测 方向 、 测 站 地 理 
位 置 、 观 测 目 标 光 谱 型 等 因素 都 影响 大 气 折射 误差 的 校正 ， 这 些 偏差 导致 几何 扭曲 模型 高 阶 项 存在 。 

几何 扭曲 模型 的 常数 项 表示 望远镜 坐标 轴 的 侦 移 程度 ， 通 常 由 于 望远镜 弯曲 随 观测 时 间 和 位 置 
发 生变 化 。 一 次 项 分 别 为 望远镜 的 比例 尺 与 取向 ， 实 际 比例 尺 的 变化 通常 由 望远镜 的 热 变化 导致 ， 模 
型 的 取向 项 用 于 对 齐 测量 坐标 系 与 标准 坐标 系 ， 与 底片 测量 仪器 的 定 标 和 高 阶 项 的 建 模 不 准确 有 关 。 
因此 ， 本 文 没有 探究 几何 扭曲 模型 的 常数 项 、 一 次 项 与 天 顶 距 的 关系 。 二 次 项 用 于 修正 底片 倾斜 导 
致 的 误差 ,通常 由 3 种 效应 引入 : 切 平面 坐标 系 之 间 的 转换 ， 从 相对 于 光 轴 倾斜 的 坐标 系 到 切 平面 坐 
标 系 的 转换 ， 以 及 光学 场 角 扭曲 (Optical Field-Angle Distortion, OFAD) 定 心 不 准确 。 三 次 项 主要 用 于 
修正 光学 场 角 扭 曲 ， 光 学 场 角 扭 曲 定 心 不 准确 会 造成 像素 坐标 位 置 偏 移 。 光 学 场 角 扭 曲 通常 表示 为 
与 坐标 有 关 的 三 次 项 形式 ， 其 微分 为 二 次 多 项 式 ， 在 光学 场 角 扭曲 系数 求解 过 程 中 产生 的 误差 在 转换 
到 切 平 面 坐标 的 过 程 中 引入 二 次 项 与 三 次 项 。 由 于 光学 场 角 扭曲 与 光线 相对 于 光 轴 的 角度 有 关 ， 我 们 
猜测 天 项 距 不 同 的 天 体 ， 其 光线 通过 望远镜 光学 系统 的 折射 角度 不 同 ， 所 以 部 分 二 次 项 和 三 次 项 的 系 
数 也 存在 与 天 顶 距 的 相关 性 。 四 次 项 可 以 修正 大 气 折射 未 完全 校正 导致 的 系统 偏 移 ， 天 顶 距 小 的 天 体 
受 大 气 折射 的 影响 小 ， 天 顶 距 大 的 天 体 受 大 气 折射 的 影响 大 ， 由 图 3 可 以 看 出 ,在 x 方向 ， 随 着 天 项 
EZK, MRH xt, cy, xy Al xy? 的 系数 对 扭曲 的 影响 逐渐 增加 ; 在 y 方 向 ， 四 次 项 x 和 wy 的 系数 
变化 同 理 。 
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几何 扭曲 模型 高 阶 项 在 * 方向 与 y 方 向 的 系数 项 呈 中 心 对 称 ， 在 * 方向， 四 次 项 xw*y 和 xy 的 系数 


与 三 次 项 x” 和 wy 的 系数 对 扭曲 的 影响 较 大 ; 在 y 方 向 ， 四 次 项 x*y 和 yy 的 系数 与 三 次 项 wy 和 Y 的 
系数 对 扭曲 的 影响 较 大 。 本 文 初步 发 现 的 规律 显示 ， 天 顶 距 越 大 ， 对 天 体 的 位 置 测量 影响 越 大 。 


4 总 结 与 展望 


本 文 探 究 了 丽江 观测 站 2. 4 m 望远镜 CCD 几何 扭曲 模型 的 变化 规律 。 我 们 选取 了 2013 年 2 月 4 
日 ~6 日 对 密集 星 场 抖动 观测 的 10 组 资料 ， 分 别 求解 几何 扭曲 模型 ， 并 用 四 阶 多 项 式 拟 合 每 个 星 
像 的 几何 扭曲 与 像素 坐标 位 置 的 关系 。 本 文 初步 发 现 几何 扭曲 模型 的 四 次 项 和 二 次 项 的 一 些 系数 与 天 
顶 距 有 一 定 的 相关 性 ， 在 * TT, WRH xi, ay, xy May, DRD Al xy 与 天 顶 距 存在 明显 的 线 
性 关系 ; 在 y 方 向 ， 四 次 项 x* 和 wy 、 三 次 项 wy 和 二 次 项 xy 与 天 顶 距 存在 明显 的 线性 关系 。 本 文 希 
望 通过 此 规律 ， 为 今后 丽江 观测 站 2. 4 m 望远镜 观测 定 标 场 节约 一 定 的 时 间 ， 通 过 较 少 次 数 的 定 标 
场 观测 推导 其 他 时 间 观 测 的 CCD 图 像 的 几何 扭曲 。 本 文 的 工作 仅 初 步 从 天 顶 距 方面 探究 几何 扭曲 模 
型 的 变化 规律 ， 在 后 续 工 作 中 ， 我 们 将 从 温度 、 气 压 等 方面 继续 深入 研究 2. 4 m 望远镜 的 几何 扭曲 变 
化 规律 。 
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Abstract: Geometric distortion is a non-negligible factor for CCD high-precision measurement, which 
usually requires careful observational work to be derived, and inevitably takes up valuable observational time. 
This paper explores the change of geometric distortion with zenith distance by using the CCD image observed 
with the GMG -2. 4 m telescope at Lijiang station. Specifically, we refer to Gaia DR3 catalog in data reduction, 
and correct the effect of differential color refraction ( DCR). Then the geometric distortion models at different 
observational times are derived and are fitted with a fourth order polynomial. The results show that the 
geometric distortion models of the CCD frames are correlated with the zenith distance ( ZD) of observational 
time. There are obvious linear relations between quadratic terms x^, x^y, x^y^, xy’, quadratic terms x^, xy 
and zenith distances in X coordinate; in Y coordinate, there are obvious linear relations between the quadratic 
terms x^, xy’, cubic term xy’, quadratic term xy of the geometric distortion models and zenith distances. 
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